Thermodynamics of the Spin-1/2 Antiferromagnetic Uniform Heisenberg Chain 
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We present a new application of the traditional thermodynamic Bethe ansatz to the spin-1/2 
antiferromagnetic uniform Heisenberg chain and derive exact nonlinear integral equations for just 
two functions describing the elementary excitations. Using this approach the magnetic susceptibility 
X and specific heat C versus temperature T are calculated to high accuracy for 5 x 10 -25 < T/ J < 5. 
The x(T) data agree very well at low T with the asymptotically exact theoretical low-T prediction 
of S. Lukyanov, Nucl. Phys. B 522, 533 (1998). The unknown coefficients of the second and third 
lowest-order logarithmic correction terms in Lukyanov's theory for C(T) are estimated from the 
C(T) data. 



The spin 5=1/2 antiferromagnetic (AF) uniform 
Heisenberg chain has a long and distinguished history 
in condensed matter physics and exhibits unusual static 
and dynamical properties unique to one-dimensional spin 
systems. It has been used as a testing ground for 
many theoretical approaches. The Hamiltonian is TL — 
JYli &i ' "S^+ii where J > is the AF Heisenberg ex- 
change interaction between nearest-neighbor spins. In 
this paper we usually set fee = 1 and gfiB = 1 where fcs 
is Boltzmann's constant, g is the spectroscopic splitting 
factor of the spins and /ib is the Bohr magneton; also, 
the reduced temperature t = T/J where T is the absolute 
temperature. 

The 5 = 1/2 Heisenberg chain is known to be exactly 
solvable H, i.e. all eigenvalues can be obtained from the 
so-called Bethe ansatz equations. Despite the amazing 
property of being integrable, the Heisenberg chain has 
defied many attempts to calculate physical observables 
including thermodynamic quantities. A rather direct 
evaluation of the partition function was constructed in 
Pj and is known as the "thermodynamic Bethe ansatz" 
(TBA), but this did not allow for high accuracy calcu- 
lations especially in the low temperature region. The 
fundamental problem in || is the necessity to deal with 
infinitely many coupled nonlinear integral equations for 
which the truncation procedures are difficult to control. 

The possibility to accurately calculate the physical 
properties of the 5 = 1/2 Heisenberg chain improved fol- 
lowing the development of the path integral formulation 
of the transfer matrix treatment of quantum systems || . 
On the basis of a Bethe ansatz solution [Q] to the quantum 
transfer matrix, Eggert, Affleck and Takahashi in 1994 
obtained numerically exact results for the magnetic sus- 
ceptibility x(t) down to much lower temperatures than 
before and compared these with their low-t results from 
conformal field theory ||. They found, remarkably, that 
x(t — ► 0) has infinite slope: their conformal field theory 
calculations showed that the leading order t dependence 



is x(t -> 0) = x(0){l + l/[21n(£ /t)]}, where the value of 
to is not predicted by the field theory. Such log terms are 
called "logarithmic corrections" . From their comparison 
of their field theory and Bethe ansatz calculations which 
extended down to t — 0.003, Eggert, Affleck and Taka- 
hashi estimated to w 7.7 M. Their numerical x(t) values 
are up to ~ 10% larger than the former Bonner-Fisher 
H extrapolation for t < 0.25. 

Lukyanov has recently presented an exact asymptotic 
field theory for x{t) and the specific heat C(t) at low 
t, including the exact value of to @- These results arc 
claimed to be exact in the sense of a renormalization 
group treatment close to a fixed point where only few op- 
erators are responsible for perturbations. Questions aris- 
ing about such calculations are whether these operators 
have been correctly identified and whether the effective 
theory has been properly evaluated. A meaningful test of 
Lukyanov's theory is only possible using numerical data 
of very high accuracy and at extremely low temperatures, 
such as we have attained in our numerical calculations to 
be presented below. 

In this Letter we present a new application of the tradi- 
tional TBA to the spin-1/2 Heisenberg chain and derive 
exact nonlinear integral equations [Eqs. (§-^) below] in- 
volving just two functions describing the elementary ex- 
citations. Our derivation evolved from earlier work by 
one of us using the powerful lattice approach By 
means of a lattice path integral representation of the fi- 
nite temperature Heisenberg chain and the formulation 
of a suitable quantum transfer matrix, a set of numeri- 
cally well-posed expressions for the free energy was de- 
rived. A serious disadvantage of this approach lies in the 
complicated and physically non-intuitive mathematical 
constructions, which strongly inhibits generalizations to 
other integrable, notably itinerant fermion models. The 
present work is a new analytic derivation of the finitely- 
many integral equations of |^|,^) by means of the intu- 
itive TBA approach. Our Eqs. (§-j|) are identical to 
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those obtained in B by a rigorous, however much more 
involved method. In our new construction, we assume 
that magnons (on paths C±) are elementary excitations 
and contain all information about the thermodynamics. 
Bound states are implicitly taken into account by use of 
the exact scattering phase probed in the analyticity strip. 
The a posteriori success of our reasoning is important for 
two reasons. First, our construction is as simple as the 
standard TBA, however avoiding the problems of dealing 
with density functions for (up to) infinitely many bound 
states. This may be of great advantage in the study of 
more complicated systems. Second, we have a simple par- 
ticle approach to the Heiscnbcrg chain which will allow 
for a study of transport properties like the Drude weight 
which has not been possible within the path integral ap- 
proach 

We also demonstrate here that using our integral equa- 
tions one can improve the accuracy and extend the tem- 
perature range of numerical calculations of %(i) and C(t) 
for the S — 1/2 Heisenberg chain on the lattice far be- 
yond those of previous calculations. We find agreement 
of our data with the above theory of Lukyanov JjJ for 
x(t) to high accuracy (< 1 x 10 -6 ) over a temperature 
range spanning 18 orders of magnitude, 5 x 10 -25 < t < 
5 x 10 -7 ; the agreement in the lower part of this temper- 
atures range is much better, 0(1O~ 7 ). For C(t), the loga- 
rithmic correction in Lukyanov's theory is insufficient to 
describe our numerical data accurately even at very low 
t, so we estimate the coefficients of the next two logarith- 
mic correction terms in his theory from our C (t) data. 
Derivation of integral equations 

We start with the partially anisotropic Hamiltonian 
H = J Ei(SfSf +1 + SfSf +1 + cos( 7 )Sf^ +1 ) - h £. Sf 
with < 7 < 7r/2 and magnetic field h. The dynamics 
of the magnons, i.e. the elementary excitations above the 
ferromagnetic state, constitute the Bethe ansatz. Mo- 
mentum p and energy e are suitably parametrized in 
terms of the spectral parameter x 



p(x) = i loj 



sinh(a; — i"//2) 
sinh(.T + ij/2) ' 



e(x) 



J — - — p (x) — n, 



(1) 



where real values are obtained for Imx = and — 7r/2, 
defining magnon bands of type "+" and "— ". 

Any two magnons with spectral parameters x and y 
scatter with phase shift Q(x — y) where 



9(z) = — ilog 



sinh(z — 27) 
sinh(z + i 7 ) 



(2) 



Next we apply the standard TBA |2| just to the magnons 
and ignore bound states! However, the magnons on 
band — are considered for spectral parameter x with 
Im x = —7 hence avoiding the branch cut in the scat- 
tering phase. 



The density functions for particles pj and holes for 
the bands j = +, — give rise to the definition of the ratio 
function rjj = p'j / Pj ■ Our analysis shows that rj + and rj- 
are analytic continuations of each other. Quantitatively 
we find r/-(x + 17) = 77+ (x) —: r)(x) subject to the non- 
linear integral equation 



1 \ £ ( x ) f 
logn{x) = — + / k 

J c 



(x-y)log(l + 7 1 - 1 (y))dy (3) 



where k{x) — w^Q'(x) and C is a contour consisting of 
the paths C+ and C_ with Im y = and —7 encircled in 
clockwise manner. Substituting log(l + f]^ 1 ) = log(l + 
j]) — log rj on the contour C+ and resolving for log 77 we 
find 



\ogr](x) = - 



e(z) 
T 



k(x - y) log(l + ii{y))dy 
R(x-y) log(l + ■n~ 1 (y))dy (4) 



with 

l{x) 

R{x) 



J—z-eoi?) 



2(tt-7) 
sinh(-| — "f)k 



2 71 " J-00 2 cosh 3rfc sinh ^-^-k 



h, e (x) = 
ikx dk 



cosh -x ' 
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(5) 



(6) 



Finally, ignoring T and h independent contributions we 
obtain the free energy as 



T 
2^ 



e (x) log [(1 + r)(x))(l +77 1 (x - 77))] dx. 



(7) 



Numerical study of low-T behavior 

Lukyanov's low-i asymptotic expansion of x(^) is 



where g(t/to) obeys the transcendental equation 
v / gexp(l/f;) = ty/t, with a unique value of to given by 
t = v /-7r/2exp(7 + 1/4) « 2.866 where 7 is Euler's con- 
stant. His expansion for the free energy per spin at h = 
yields the specific heat per spin as 



at, g (t) = § 



l + ^ 3 + 0(.9 4 ) 



(2)3 5 / 2 t 3 



57T 



[1 + 0{g)] , 



(9) 



where the exact prefactor 2t/3 was found by Affleck in 
1986 [[To), and the prefactor 3/8 in the logarithmic cor- 
rection term agrees with p|0,0 . 
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Numerical data for %(*) and C(t) were obtained using 
our free energy expression (|7|). These data are consid- 
erably more accurate than those presented previously in 
[|L Our x(t)J data, and the exact value l/n 2 at t = 
[ 13| , are plotted in Fig. [l| The calculations have an ab- 
solute accuracy of sw 1 x 10~ 9 . The data show a maxi- 
mum at a temperature * max = 0.6 408 510(4) with a value 
X max J = 0.146 926 279(1), yielding the J-independent 
product x max T max = 0.0 941579(1). These values are 
consistent within the errors with those found by Eggert, 
Affleck and Takahashi ||, but are much more accurate. 

The differences between our low-* Bethe ansatz %(<) J 
calculations and Lukyanov's theoretical xit,g(*)<^ predic- 
tion in Eq. (Q) are shown in Fig. |^. The error bar on each 
data point is the estimated uncertainty in Xit,gJ arising 
from the presence of the unknown 0(g i ) and higher-order 
terms in Eq. (||), which was arbitrarily set to g 4 (t)/ir 2 ; 
the uncertainty in the t 2 contribution, ~ ^/3t 2 g(t)/-K 3 , is 
negligible at low * compared to this. At the lower tem- 
peratures, the data agree extremely well with the predic- 
tion of Lukyanov's theory. At the highest temperatures, 
higher order t n terms also become important. Irrespec- 
tive of these uncertainties in the theoretical prediction at 
high temperatures, we can safely conclude directly from 
Fig. ^ that our numerical x(t) data are in agreement 
with the theory of Lukyanov to within an absolute 
accuracy of 1 x 10~ 6 (relative accuracy 10 ppm) from 
* = 5 x 10~ 25 to t = 5 x 10~ 7 . The agreement at the 
lower temperatures, 0(1O -7 ), is much better than this. 

Our C(t) data for t < 2 are shown in the in- 
set of Fig. U and have an estimated accuracy of 3 x 
10 _10 C(*). The data show a maximum with a value 
C max = 0.3 497121235(2) at a temperature *£ ax = 
0.48 028 487(1). The electronic specific heat coefficient 
C(t)/t is plotted in Fig. || These data exhibit a max- 
imum with a value (C/*) max = 0.8 973 651576(5) at 
i'° ax = 0.30 716 996(2). The existence of low-* log cor- 
rections to C (*) is revealed in the top plot of AC (*) /* in 
Fig. |, where AC(t) = C(t) - 2*/3 and 2*/3 is the low-* 
limit of C(t). The influence of the g 3 log correction term 
in Eq. (||) is evaluated by subtracting it in the plot of 
AC(t)/t as shown by the middle curve in Fig. |. The 
t = singularity is still present but with reduced ampli- 
tude; this demonstrates that additional logarithmic cor- 
rection terms are important within the accuracy of the 
data. 

We estimate the unknown coefficients of the next two 
logarithmic correction (g 4 , g 5 ) terms in Eq. (||) from our 
C(t) data as follows. From Eq. (||), if we plot the data 
as [C{t)/t- (2/3)(l + 3.g 3 /8)]/5 4 vs g and fit the lowest-* 
data by a straight line, the y-intercept gives the coeffi- 
cient of the g A term and the slope gives the coefficient 
of the g 5 term. We fitted a straight line to the data in 
such a plot for 5 x 10~ 25 < * < 5 x 10~ 9 as shown by the 
weighted linear fit in Fig. where the parameters of the 



fit are given in the figure. By subtracting the influences 
of these two logarithmic correction terms from the mid- 
dle data set as shown in the bottom data set in Fig. |J, 
the singular behavior as * — > is largely removed, leaving 
a behavior which is close to a t 2 dependence as predicted 
by the last term in Eq. (||). Further discussion of the 
predictions of H, and high-accuracy fits (0 < * < 5) to 
our C(*) and x(*) data and the respective exact * = 
values, will be presented elsewhere Q. 

In conclusion, we have presented an analytic approach 
to the thermodynamics of the 5 = 1/2 AF Heisenberg 
chain on the basis of a finite number of elementary excita- 
tions. We envisage that this approach can be generalized 
to study a variety of other systems such as Hubbard and 
t-J models, quantum spin chains with higher symmetries 
and systems with orbital degrees of freedom. Our free 
energy expression has allowed numerical calculations of 
x(t) and C(f) for the Heisenberg chain to be carried out 
to much higher accuracy and to much lower temperatures 
than heretofore attained. Our %(*) data are in excellent 
agreement with the theory of Lukyanov |?J at low *. The 
logarithmic correction in Lukyanov's theory for C{t) is 
found insufficient to describe our C(*) data accurately 
even at very low *. However, the * dependence of the de- 
viation agrees with the form of his theory, which enabled 
us to estimate the unknown coefficients of the next two 
logarithmic correction terms in his theory for C(t) from 
our C(t) data. Thus we have verified Lukyanov's theory 
[Q of a critical system perturbed by marginal operators 
and have given evidence that his asymptotic expansion 
can be systematically extended to higher order. 
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FIG. 3. Electronic specific heat coefficient C/T versus tem- 
perature T for the 5* = 1/2 AF uniform Heisenberg chain. In 
the inset the specific heat C versus T is shown. 
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FIG. 1. Magnetic susceptibility \ at low temperature T for 
the spin 5=1/2 antiferromagnetic uniform Heisenberg chain. 
In the inset x(^) is shown on a larger temperature scale. 
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FIG. 2. Semilog plot vs temperature t at low t of the dif- 
ference between our Bethe ansatz magnetic susceptibility \J 
data and the prediction Xit.g^ of Lukyanov's theory Q. The 
error bars are the estimated uncertainties in X\t,g{t)J ■ 
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FIG. 4. Difference AC/t between our Bethe ansatz elec- 
tronic specific heat coefficient data and the exact coefficient 
2/3 at t = (top data set), versus temperature t. Succes- 
sive data sets show the influence of subtracting cumulative 
logarithmic correction terms. 



Fit: 1.025(1) + 2.08(2) g 
Fit range: 0.017 < g < 0.046 
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FIG. 5. Plot showing the estimation of the coefficients of 
the g 4 and g 5 logarithmic correction terms in Eq. (|^). The 
error bars on the data points are smaller than the symbols. 
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